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Abstract 

The large decay rate observed by Belle for B + — > -0(377O)-RT + , which is comparable 
<^ to B+ -» V(3686)iT+, 

might indicate either an unexpectedly large S-D mixing angle 
' |0| « 40° or the leading role of the color-octet mechanism in D-wave charmonium 

production in B decay. By calculating the production rate of "0(3770) in the continuum 
\ e + e~ annihilation at y/s = 10.6 GeV with these two possible approaches (i.e. the large 

' S-D mixing and the color-octet mechanism), we show that the measurement for this 

process at Belle and BaBar may provide a clear cut clarification for the two approaches. 
In addition, the radiative El transition ratio T(-0(377O) — > 7Xc2)/T (-0(3770) — ► 7Xcl) 
I/"") | may dramatically change from ~ 0.04 (for 6 ~ 0°) to ~ 200 (for 6 ~ —40°) due to the 

large S-D interference effect, thus the El transition measurement of -0(3770) at BES 
. and CLEO-c will also be very useful in clarifying this issue. 

PACS number(s): 13.66.Bc, 14.40. Gx, 12.38.Bx 

a: 

1 Introduction 

>- ! The S-D mixing for ip' = -0(3686) and -0" = -0(3770) is of great interest in charmonium 
physics. If we neglect the charmed meson pair component which is due to coupling to decay 
?-h ' channels, -0' and -0" may be approximately expressed as 

lip'} =cos6\2 3 S 1 } + stn6\l 3 D 1 }, \ ip") = cos9 | l 3 £>i) + sind \2 3 S 1 ). (1) 

A rough estimate of the S-D mixing angle may be obtained by using the ratio of the observed 
leptonic decay widths pQ of -0' and ip" and neglecting the D-wave component contribution to 
these leptonic decay widths tan 2 # = ffef^pp) ~ 0.12, which results in 9 w ±19°. However, 
if the D-wave contribution to leptonic decay widths is further included, potential model 
calculations e.g. in IU give two solutions: 9 w -10° to -13° or 9 « +30° to +26° 
(6 w -10° and +30° in 0, 9 w -13° and +26° in and 9 w -12° and +27° in 0]). The 
small mixing solution (i.e. 9 ~ —10°) is compatible with the results obtained in models with 
coupled decay channels jSHH]- Moreover, 9 ~ —10° is favored by the ip' — > ^Xcj data whereas 
9 ps +30° would lead to T(ip' — > 7X c o) ~135 KeV due to a large positive S-D interference [3], 
which is higher than its observed value by a factor of 6 and is therefore disfavored. The S-D 
mixing may have many interesting phenomenological consequences. It might slightly |2] or 
substantially [7||H] affect the -0(3770) — > J/tpnir decay , which has been recently observed by 
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BES [9 . It would also have effects on the h c search via the ip' — > h c 7c° decay ^U]- The S-D 
mixing could even provide an explanation for the notorious pn puzzle that the suppression 
of -0' — ► pn is due to a destructive interference between the S and D wave states jlj, and it 
might also be useful in explaining JI] the recent observed enhancement of -0' — > K^Ks by 
BES [12]. Moreover, the ratio of jl^'j may sensitively change from 0.04 (for 9 « 0°) 

to 0.22 (for 9 « —10°) and to 0.06 (for 9 « +30°) 0). The experimental examination of 
these radiative transitions for ^(3770) at BES and CLEO-c in the near future will be an 
interesting test for the S-D mixing. 

Recently Belle Collaboration ^3] has observed -0(3770) for the first time in the B meson 
decay B + -> ^(3770)^+ with a branching ratio of (0.48 ± 0.11 ± 0.07) x 10~ 3 , which is 
comparable to B(B + — > i/;'(368d)K + ) = (0.66 ± 0.06) x 10" 3 [T]. This is quite surprising, 
since conventionally 0(3770) and ■0(3686) are regarded as mainly the l 3 Di and 2 3 1 S 1 color- 
singlet cc states respectively, and the coupling of l 3 Di to the cc vector current in the weak 
decay effective hamiltonian is much weaker than that of 2 3 S\. One possible explanation is 
that the S-D mixing for -0(3770) and 0(3686) is very large, much larger than previously 
expected . If in the B meson decay we neglect the l 3 Di contribution, which is expected to 
be much smaller than the 2 3 S' 1 contribution, we would get the S-D mixing angle \9\ « 40.4° 
from the observed decay rate ratio 

g(g + ^ V>(3770)Ar+) _ 2 

B{B+ -> ^(3686)K+) " ^ °~ °' 73 - (2) 

Although this large mixing angle seems to be not compatible with all our previous knowledge 
about the S-D mixing, it is still worthwhile to test it with new experiments. Another possible 
explanation is that the color-octet mechanism in nonrelativistic QCD (NRQCD) may play 
the leading role in the D-wave charmonium production in B meson decays, and it was 
predicted jHj that for a pure D-wave state the branching ratio B{B — > -0(377O)X) rs 0.28%, 
which is comparable to £>(£> — > -0(3686)X) = (0.35 ± 0.05)% [T] (see also [T^j for similar 
discussions). In a series of papers [TBJ H3 EH EE] it was pointed out that based on the Fock 
state expansion and velocity scaling rules of NRQCD ^H] the production rates of D-wave 
charmonium states including the 3 Di 0(3770), which would predominantly decay to the DD 
meson pair, and an expected narrow 3 D2 state, which could have some decay fraction to 
J/ip7i + 7r~, would be comparable to that of J/ip and if)(2S) in the Z° decay [16J, the pp 
collision at the Tevatron ^7], the B meson decay [Tl], and the fixed target experiments 
18 , whereas in the conventional color-singlet model the D-wave production rates, which are 
proportional to the squared second derivative of the cc wave function at the origin, should 
be greatly suppressed. Despite of certain uncertainties associated with the values of color- 
octet matrix elements, the color-octet contributions are expected to be dominant, which are 
larger than the color-singlet contributions by more than one order of magnitude in those 
processes. For instance, in the large pt charmonium production at the Tevatron, the 3 Dj 
states could have large rates, comparable to i/;(2S), because in both cases the dominant 
production mechanism is expected to be gluon fragmentation into the color-octet 3 S± cc 
intermediate state which then evolves respectively into the physical color-singlet 3 Dj and 
ip(2S) states by emitting two soft gluons via double El transitions with the same order 
transition probabilities |17j . Similarly, the -0(3770) production in the B meson decay could 
provide another interesting test for the color-octet mechanism in NRQCD. 
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In order to distinguish between the large S-D mixing and the NRQCD color-octet mech- 
anism in the ?/>(3770) production in the B meson decay, we suggest measuring the -0(3770) 
production in continuum e + e~ annihilation at the Belle and BaBar energy ^/s = 10.6 GeV, 
and we will give calculations of the production cross sections in the following sections. We 
will show that the ■0(3770) production in e + e~ annihilation via double cc does not receive 
large contributions from the color-octet mechanism but is very sensitive to the S-D mixing. 

J/ip inclusive production in e + e~ annihilation has been investigated within the color- 
singlet model and the color-octet model |22 • BaBar [21] and Belle j2S] have 
measured J/ip production rate in continuum e + e~ annihilations at ^/s = 10.6GeV^, which is 
found to be much larger than the color-singlet prediction. More interestingly, Belle further 
finds J/ip production to be dominated by the double cc production j2H]- The measured 
exclusive cross section for e + + e~ — > J/ip + 7] c is an order of magnitude larger than the 
theoretical values[2Zj, and the measured inclusive cross section for e + + e~ — > J/tp + c + c is 
more than five times larger than NRQCD predictions. Therefore the inclusive charmonium 
production via double cc is particularly interesting and worth investigating. As in the case 
of inclusive production of J/ip, rj c , and Xcj(J=0, 1, 2)[2E], here we will also concentrate on 
the double cc production for the D-wave charmonium states 5j( J = 1, 2, 3). 



2 Color-singlet contribution to 5\ ( 3 Di) charmonium 
production via double cc in e + e~ annihilation 

In NRQCD the Fock state expansion for the D-wave (without S-D mixing) charmonium 
5j(J= 1,2,3) is 

\Sj) =0(l)\cc( 3 Dj,l)) 
+0(v)\ccCPj,,8)g) 

+0(v 2 )\cc( 3 S 1 ,8) gg) + ---. (3) 

Following the NRQCD factorization formalism, the scattering amplitude of the process 
e (Pi) + e+ G°2) - ► 7* — ► cc( 25+1 Lj 1 ' 8a ' ) )(p) + c(p c ) + c(p 5 ) in Fig. [T] is given by 

A(e-( Pl ) + e + (p 2 ) - ccC S+1 L { j M) )(p) + c{p c ) + c(p 5 )) = ^EEE 

L Z S Z «i«2 jk 

x (si,s 2 | SS Z )(LL Z ; SS Z | JJ z )(3j;3k \ 1,8a) 

f A(e-( Pl ) + e + (p 2 ) -> Cjif, s x ) + c fc (f ; s 2 ) + c,(f ; s 3 ) + c 4 (f ; s 4 )) (L = S), 

x e ;(L z )^(e-(p 1 ) + e+(p 2 )^c J (f; Sl ) + c fc (f; S2 ) + Q(f;s 3 ) + Q(f; S 4)) (L = P), 

{ \e* ap {L z )A^{e-{ Pl ) + e+(p 2 ) - Cj -(f ; Sl ) + c fe (|; s 2 ) + Q (f ; s 3 ) + P.(f ; s*)) (i = D). 

(4) 

where cc( 25+1 L?' 8a) ) is the cc pair produced at short distances, which subsequently evolve 
into a specific charmonium state at long distances, A a and A a @ are the derivatives of the 
amplitude with respect to the relative momentum between the quark and anti-quark in the 
bound state. For the case of color-singlet state, the coefficient Cl can be related to the origin 
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of the radial wave function (or its derivatives) of the bound state as 



1 ^ 

Cd = £ I i£(0) | 2 • (5) 

The spin projection operators and their derivatives with respect to the relative momentum 
are 

P 1Sz (p,0) = -L 4(S z )(^ + 2m c ), (6) 
P is>'°) = ii^ a f{S z ){t> + 2m c )-{t>-2m c ) 4(S,)>f\. (J) 

<(p,o) = ^^h a ns*)-f + -f Ks z )i a ]. (8) 

The calculation of cross sections for e~ + e + — > 7* — > charmonium +cc is straightforward. 
As in Ref. |2D1 we write the differential cross section as 

da(e + + e~ — > 7* — > charmonium + cc) 4Ci)Q; 2 af , , . a(z) . , . 

1 = — qi W{z) H —J. (9) 

az 81m c 3 

where z = 2E^/y/s. The expressions of S(z) and a(z) for 5j are lengthy and will be given 
in the appendix for J = 1,2. With Eq. Q we can evaluate the inclusive cross sections for 
8 j. The input parameters used in the numerical calculations are 

m c = 1.5GeV, a s (2m c ) = 0.26, a = 1/137, | R%(0) \ 2 = 0.015Ge1/ 7 |2ni. (10) 

and the obtained cross section for the 5i at ^/s = 10.6 GeV is 

a(e + + e~ -> 7* -> 5i + cc) = 2.5 fb. (11) 
Here we also give the calculated cross section for the S 2 at y/s = 10.6 GeV 

a(e + + e~ -> 7* -> 5 2 + cc) = 2.4 fb. (12) 



We also find that as in the case of other charmonium states j28] the calculated cross sections 
for 5\ in eq. (jllj) and 62 in eq. ()12j) are substantially smaller than those obtained in the 
fragmentation approximation at y/s = 10.6 GeV which would cause a enhancement factor 
of 1.5 and 2.3 respectively. 

Comparing eq.(JTTJ) with a(e + + e~ — > 7* — > J /ip + cc) = 148 fb calculated for the J /ip 
in [2H], we see that the inclusive double cc cross section for the D-wave 1 state is smaller 
than that for the S-wave 1 states by a factor of 60. This illustrates the expectation that 
within the color-singlet model the suppression of D-wave state production relative to the 
S-wave state production is usually about two orders of magnitude. 
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3 Color-octet contribution to Si production via double 
cc in e + e~ annihilation 

The color-octet contributions to 8\ production via double cc in e + e~ annihilation come from 
the Feynman diagrams in Fig ^ and Fig El According to the NRQCD factorization and 
velocity scaling rules^Hl, the contributions of the second term and third term in the Fock 
state expansion of 8\ in eq. (jSJ) are of the same order in the quark relative velocity v as the 
corresponding terms in the Fock state expansion of J/ip or ip(2S). As a rough estimate for 
the nonperturbative matrix elements we may choose 



Here the color-octet matrix elements for J/tp were extracted from the J/ip data at the 
Tevatron (see Ref. [SHI Ell f° r detailed discussions). There are large uncertainties with < 
Og^( 3 P ) > /m 2 c and < Og^( 1 S , ) > /3 and their combinations, and here we have assumed 
that they are equal and take the largest fitted values from jHSj to avoid underestimates of 
the color-octet contributions, (note that these two matrix elements may be overestimated 



In Fig. Q the color-octet contribution can be obtained from the corresponding color- 



from the color-octet 3 S U 3 P 3 P x 3 P 2 states to be 0.12fb, l.lfb, 0.29fb, 0.14fb respectively. 
The total color-octet contribution to the cross section from Fig. ^ is 1.65fb. 

In Fig. |21 the color-octet contributions come from four different (the upper two and the 
lower two) diagrams. The upper diagrams only contribute via the color-octet 3 Si state, and 
the differential cross section reads 



< OfreSt) >^< C 8 J/ ^( 3 5i) >= 1.06 x 1(T 2 GeV 3 , 

< O^CPq) > /ml ^< O J 8 /lp ( 3 P ) > /m 2 c = l.lx 10" 2 GeV 3 , 

< OfrCSo) >« O^CSo) >= 3.3 x 1(T 2 GeV 3 , 

< Of ( 3 Pj) >= (2J + 1) < Og( 3 P ) > . 



(13) 
(14) 
(15) 
(16) 




da. 



octet 



16a 2 a 2 < O s 8 1 ( 3 S 1 ) > 



M 



2 



(17) 



dz 



27 m c 



where I M I 2 takes the form 




[3S 4 - I25 2 (z - 2) + 16(10 + z{z - 10))] + 
(z - 2) 2 [35 4 - 85 2 (3z - 4) + 32(2 + z(z - 2))] 




(18) 
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The numerical result for the cross section is 

c?octet(e + e- -> 5 lC c) = 4.5 fb. (19) 

The lower diagrams in Fig. |2] make contributions to 5% production via the color-octet 
3 Pj( J = 0, 1, 2) and 1 5'o intermediate states. We find the following results. 



da octet 32a 2 a 2 < 0* ( 25+1 A/) > . /(z 2 - 5 2 )(1 - z) , -,. 2S+1t . |2 



dz 27m r V 4 + <5 2 - Az 



M( 2b+1 Lj) \\ (20) 



- 3 2 _ (16tt(8 + 3£ 2 - 8z)(S 6 + 8£ 2 (-4 + z)z 2 + 8^ 4 + 2£ 4 (18 + (-10 + z)z))) 
1 0)1 (275 2 s 3 (A + 5 2 -Az) 2 (-2 + z y) 1 ' 



M( 3 P] 



3n \ |2 



(95 2 s 3 (4 + 5 2 - 4z) 2 (-2 + z)*) 
x (16tt(8 + 35 2 - 8z)(5 6 + 16(-2 + 2) V + 2<5 4 (18 + (-12 + z)z) 
- 4<5 2 (-2 + z)(16 + z(-20 + 3z)))) (22) 



- 3 2 _ (32tt(8 + 35 2 - 8z)(5* + 2£ 4 (-3 + z)(-2 + z) + 4^ 4 - 2b 2 z\2 + z))) 
1 2j 1 (95 2 s 3 (4 + <5 2 -4z) 2 (-2 + ^) 4 ) 1 ' 



I w r i „ x ,2_ (-87r(8 + 3J 2 -8z)(J -*)(* + *)) 

1 M( 5o) 1 " (9s 2 (4 + 5 2 -4^) 2 (-2 + ^) 2 ) (24) 

With the chosen parameters mentioned above, the contributions to the 8\cc production 
cross section come from the color-octet 3 Pq, 3 Pi and 3 P2 are 0.18fb, 2.7fb and 0.87fb respec- 
tively. The contribution of the color-octet 1 So is negligible. 

The total color-octet contribution to the cross section from Fig. |2]is 8.3fb, and the sum 
of the color-octet contributions from Fig. ^ and Fig. |2]is 9.9fb. This number could become 
substantially smaller if we use the matrix element values given in |34j . In any case, the color- 
octet contribution to the D-wave states production cross sections should be of the same order 
as the color-singlet contribution, because from the 5j Fock state expansion in eq. and from 
Fig. Q and Fig. |2] it is easy to see that the color-octet and color-singlet contributions are of 
the same order in both short distance (0(a 2 )) and long distance (0(v 4 )) parts. Based on the 
above calculations including both the color-singlet and color-octet contributions, we think 

or(e + + e~ -> 7* -> <5i + cc) « 10 fb (25) 

should be a reasonable estimate in NRQCD for the ^(3770) production rate when the S-D 
mixing is neglected. 
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4 S-D mixing and ^(3770) production in e + e annihila- 
tion and B decay 

The calculated 0>(377O)cc production rate could be significantly enhanced by the S-D mixing 
if the mixing angle is large. With the calculated J/ip production cross section a(e + + e~ — > 
7* — > J/ip + cc) = 148 fb |2H], the ip(2S) production cross section can be approximately 
obtained by the scale factor | i?2s(0) | 2 / | -Ris(O) | 2 , and then we have a(e + + e~ — > 7* — ► 
ip{2S) + cc) = 90 /fo (the color-octet contributions to the S-wave charmonia are small and 
negligible). If this estimate for the ip(2S) makes sense, with the large S-D mixing angle 
9 pa ±40° we would have 

a(e + + e~ -> 0(3770) + cc) pa a{e + + e~ -> ^(2£) + cc) x tan 2 « 66 fb. (26) 

This value is much larger than 10 fb, the value in eq.(|25p obtained without S-D mixing in 
NRQCD. 

However, as we already mentioned, the Belle observed double charm production cross 
section for the J/ip has a much higher value [2l>] 

a(e + + e~ -> J/ip + cc) = (0.87 ± 0.15 ± 0.12) pb, (27) 

which is larger than the theoretical expectation in NRQCD by a factor of 6 (as a rough 
estimate we neglect the feed down contribution from the ip(2S) and \cj states in the J/ip 
production cross section). If this also happens for the ip(2S), which is very likely, we would 
expect 

a(e + + e" -> ip(2S) + cc) pa 530 fb (28) 

to be the observed double charm production cross section for the ip(2S). Then with the 
large S-D mixing angle 9 pa ±40° we would have 

a(e + + e~ -> ^(3770) + cc) pa a(e + + e~ -> ^(25) + cc) x tan 2 w 387 fb. (29) 

This is more than an order of magnitude larger than the value in eq. ()25)1 obtained without 
S-D mixing in NRQCD. 

In any case, we see that the large S-D mixing angle 9 pa ±40° would result in a much 
higher cross section for the ^(3770) double charm production than that without S-D mixing 
or with small S-D mixing angle like 9 ~ —10°. This is due to the fact that the double charm 
production rate for a pure 2 3 S\ state is much higher than that for a pure l 3 Di state in 
NRQCD. We therefore suggest measuring the production cross sections of -0(3770) + cc and 
0(3686) + cc in the continuum e + e~ annihilation at ^Js = 10.6 GeV by Belle and BaBar. 

In contrast, in the B inclusive decay to 0(3770), the decay rate is expected to be insensi- 
tive to the S-D mixing in NRQCD, because both the 2 3 S± and l 3 Di final states produced in 
B decay are dominated by the color-octet intermediate states with the same quantum num- 
bers in the two cases and thus they may have comparable production rates. The color-octet 
dominance in B decay relies on two observations. The first is that in the effective weak in- 
teractions the squared short distance coefficient at the b quark mass scale for the color-octet 
part is larger than that for the color-singlet part by more than an order of magnitude. The 
second is that the color-singlet S-wave contribution is further suppressed by QCD radiative 
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corrections and therefore negligible (see [US] [US] [SZj for detailed discussions). Therefore, de- 
spite of certain uncertainties related to the values of color-octet matrix elements and other 
parameters, the qualitative features for the D-wave charmonium production in e + e~ annihi- 
lation and B decay should hold and be tested. 



5 S-D mixing and ^(3770) El transitions 

As mentioned already, the small S-D mixing angle like 9 ~ —10° is favored by the observed 
El transition rates of "0(3686). However, as another independent check for the mixing angle, 
measurements on the ■0(3770) El transitions will be also very useful. In ref. 3\ the El 
transition rates of ■0(3770) were calculated for 9 = 0°, —10°, +30°, and the S-D interference 
effects were found to be significant. Based on the same potential model we now estimate the 
El transition widths for 9 = ±40.4° and find 

T(^(3770) -> 7% c j) = 386, 0.32, 66 KeV (30) 

for 3=0,1,2 with 9 = -40.4°; and 

T(^(3770) -> 7x c j) = 52, 203, 28 KeV (31) 

for J=0,l,2 with 9 = +40.4°. 

We see that the S-D interference effects in the -0(3770) El transitions are essential. In 
particular, the ratio 

r(0(377O) -> 7Xc 2 ) 
1^(3770) - 7X ci) 

will dramatically change from ~ 0.04 (for 9 pa 0°) to ~ 200 (for 9 ~ —40°) due to the large 
S-D interference effect. Similar discussions for 9 pa —12° and +27° can also be found in j^j. 

We hope these measurements can be performed at BES and CLEO-c in the near future. 
They will be very helpful in clarifying the S-D mixing problem for 0(3770) and 0(3686). 



R 2/1 = ^):)ZZl i (32) 



6 Discussions and summary 

We first discuss about the uncertainties associated with the color-octet matrix elements with 
different choices from eqs.(13)-(15). According to the NRQCD velocity scaling rules, we may 
have 

(Ol'CDj)) ~ m c V, (0 8 5j ( 3 Pi)> ~ my, (Ot'CS^) ~ my. (33) 

If we use (Ct( 3 Pi)) as the input parameter, (O^^Dx)) = 4gfe | P£(0) \ 2 = O.lQGeV 7 , we 
would have the following matrix elements 

tfW B |fl?W B |fl?W> = om (34) 

m' c m\ ml 

Then the color-octet contributions to the cross section cx(e + e~ — > ( 2S+1 Lj) + cc — > 8% + cc) 
can be estimated to be 0.36 fb, 1.1 fb, 0.28 fb, 0.12 fb for 3 S{ 8) , 3 P (8) , 3 P^, 3 P 2 (8) color-octet 
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intermediate states respectively from Fig. 1; and 13.5 fb, 0.17 fb, 2.6 fb, 0.83 fb respectively 
from Fig. 2. The total contribution to the cross section cr(e + e - — > 7* — > 5\ + cc) is about 
20 fb, a factor of 2 larger than 10 fb, the value given in eq.(25) obtained by using eqs.(13)- 
(15). The above results may demonstrate the possible uncertainties associated with the 
color-octet matrix elements, which are estimated according to the NRQCD velocity scaling 
rules but with different choices for the input parameters. We expect that these uncertainties 
do not change our main analysis that the -0(3770) production in e + e~ annihilation receives 
much smaller contributions from the color-octet channels than from the large S-D mixing. 

We suggest performing the measurements on the 4>(2S) and -0(3770) production in the 
continuum e + e~ annihilation by Belle and BaBar in order to check (i) if the ip(2S) +cc cross 
section is as large as about 0.5pb, which would confirm the S-wave (not only for IS but 
also for 2S) charmonium production enhancement via double cc; (ii) if the cross section of 
-0(3770) + cc is comparable to that of ip(2S) + cc, which would confirm the large S-D mixing 
angle; (iii) if the cross section of -0(3770) + cc is as small as say 10 fb, which would favor 
the prediction in NRQCD. It is very likely that Belle and BaBar will find a strong signal 
for 0(2,5) + cc, because experimentally Belle has found a large cross section for the inclusive 
ip(2S) +X cross section, which is comparable to that for J/ip + X j2H]- If a large ip(2S) +cc 
cross section turned to be the case, then one can easily distinguish between the large S-D 
mixing and the NRQCD prediction (together with the small S-D mixing) by simultaneously 
measuring the -0(3770) -f cc cross section. 

In conclusion, we notice that the large decay rate observed by Belle for B + — ► -0(3770) K + , 
which is comparable to B + — > tp(3686)K + , might indicate either an unexpectedly large 
S-D mixing angle \9\ ~ 40° or the leading role of the color-octet mechanism in D-wave 
charmonium production in B decay. By calculating the production rate of 0(3770) in the 
continuum e + e~ annihilation at \/s = 10.6 GeV with these two possible approaches (i.e. the 
large S-D mixing and the color-octet mechanism), we show that the measurement for this 
process at Belle and BaBar may provide a clear cut clarification for the two approaches. 
In addition, the radiative El transition ratio r(0(377O) — > 7Xc2)/r(0(377O) — > jXci) ma J 
dramatically change from ~ 0.04 (for 6 ps 0°) to ~ 200 (for 6 rs —40°) due to the large S-D 
interference effect, thus the El transition measurement of ^(3770) at BES and CLEO-c will 
also be very useful in clarifying this issue. 
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Appendix 

In this appendix, for 5i and 5 2 we give the expressions of S(z) and a(z) which are defined 
in Eq. Q. 



x 



X 



q 16* (l-z){z*-P) 

Sl "" 2255 6 s*(5-z)(-2 + zy o z 7 (5 + zy ^ A + 5 2 - Az 



-55296005 12 - 13824005 14 + (423936005 10 + 382464005 12 + 55296005 14 ) 



(1246003205 8 + 2602598405 10 + 1092556805 12 + 101145605 14 U 2 
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+ (1459814405 6 + 7028736005 s + 6870374405 10 + 1754995205 12 

+ 109555205 14 ),2 3 - (1061683205 4 + 7329792005 6 + 17068032005 s 

+ 10373913605 10 + 1788595205 12 + 77414405 14 ),2 4 + (4836556805 4 

+ 15027609605 6 + 23246745605 s + 9920908805 10 + 1216185605 12 

+ 36720005 14 ),z 5 - (1061683205 2 + 7591526405 4 + 14484213765 6 

+ 18669870085 s + 6066188805 10 + 545281925 12 + 11101445 14 )z 6 

+ (5249433605 2 + 985907205 4 + 3293163525 6 + 8035307525 8 + 2217473285 

+ 152457605 12 + 1886885 14 )/ + (70778880 - 14362214405 2 + 10350428165 4 

+ 4486123525 6 - 1463897605 s - 420819205 10 - 20670245 12 + 266165 14 )/ 

- (566231040 - 27428454405 2 + 10700042245 4 + 823208965 6 - 403855365 

- 732805 10 + 2629445 12 + 31805 14 )^ 9 + (2235432960 - 31455805445 
+ 4930457605 4 - 2871147525 6 - 344988805 s + 18112805 10 - 411605 
+ 6305 14 )^ 10 - (5315624960 - 11540561925 2 + 7735485445 4 - 1013722885 

- 120394245 s - 2589405 10 + 54455 12 )z n + (8075444224 + 16566722565 2 
+ 11232215045 4 + 503336645 6 - 11964965 s + 1108625 10 - 3155 12 )^ 12 

- (8028487680 + 24105758725 2 + 7237817605 4 + 381822565 6 + 7919605 s 

- 7955 10 )^ 13 + (5278310400 + 13897953285 2 + 2365640325 4 + 85830485 6 
+ 163025 8 )z 14 - (2310758400 + 4454292485 2 + 398742405 4 + 5909845 6 )z 
+ (673812480 + 846069765 2 + 35189765 4 )z 16 - (124616704 + 90823685 2 )z 

, z^A + 5 2 - Az + 2 J (I - z){z 2 - 5 2 ) 

+ 11939840z 18 ] + 155 2 (z-2) 4 ln r 



rl2 



15 
17 



zy/A + 5 2 -Az- 2^(1 - z){z 2 - 5 2 
[230405 12 - (1766405 10 + 691205 12 )^ + (5191685 s + 5859845 10 



+ 898565 1 2)z 2 - (6082565 6 + 16204805 s + 7992965 10 + 639365 12 



(4423685 4 + 13808645 6 + 19129605 s + 5789445 10 + 268805 m ~ 4 



z 



3 



12^6 



- (7249925 4 + 3781125 6 + 9244805 s + 2272485 10 + 61445 12 )z 5 

- (10321925 2 + 2560005 4 + 14376965 6 + 221765 s - 469765 10 - 10565 
+ (23101445 2 + 5038085 4 + 12625285 6 + 1088165 s - 119165 10 - 2525 12 )^ 7 

- (2293765 2 - 8565765 4 - 1949445 6 - 299765 s - 35925 10 - 425 12 )^ 8 

+ (294912 - 34283525 2 - 13972485 4 - 5711845 6 - 255165 s - 4475 10 )^ 9 



10 



- (1392640 - 42598405 2 - 5836805 4 - 2150725 6 - 30785 s + 215 l0 )z 10 
+ (1998848 - 27207685 2 + 207365 4 - 160525 6 + 3455 8 )^ n 

- (1163264 - 10319365 2 + 988165 4 + 35545 6 )^ 12 + (284672 - 1728005 2 
+ 222885 4 )^ 13 - (35840 - 53765 2 )^ 14 + 4096^ 15 ]}. 



+ 



16tt I (1 - z)(z 2 - S 2 ) 

2255 6 s 4 (5 - z){-2 + zY°z 7 {5 + z) *~ Z \ A + 5 2 -Az 
x [55296005 12 + 13824005 14 - (460800005 10 + 391680005 12 + 55296005 14 )^ 
(1246003205 s + 2488320005 10 + 1054771205 12 + 101145605 14 )^ 2 
(486604805 6 + 6131712005 s + 5746329605 10 + 1552627205 12 
+ 109555205 14 )z 3 + (353894405 4 + 780288005 6 + 13091635205 s 
+ 7554892805 10 + 1409510405 12 + 77414405 14 )^ 4 - (1179648005 4 

- 3594854405 6 + 15671808005 s + 6205248005 10 + 821395205 12 
+ 36720005 14 )z 5 + (3185049605 2 + 2907340805 4 - 15875706885 6 
+ 10710246405 s + 3243015685 10 + 295224965 12 + 11101445 14 )/ 

- (10321920005 2 + 2356428805 4 - 31950499845 6 + 2564321285 s 

+ 974442245 10 + 47994245 12 + 1886885 14 )/ - (-70778880 + 2762342405 2 
+ 13703577605 4 + 43140567045 6 + 2825500165 s - 66699525 10 + 6473765 12 
+ 266165 14 )/ + (-566231040 + 65295155205 2 + 47111557125 4 
+ 41887114245 6 + 3555712005 s + 62149765 10 + 8323205 12 + 31805 14 )/ 

- (-2235432960 + 143127183365 2 + 68292997125 4 + 28443735045 6 

+ 1998980485 s + 41644325 10 + 426965 12 + 6305 14 )^ 10 + (-5315624960 
+ 159987466245 2 + 55888092165 4 + 13295578885 6 + 749368965 s 
2659005 10 + 195155 12 )^ n - (-8075444224 + 106532003845 2 
28559093765 4 + 4400174405 6 + 162533285 s + 2848505 10 + 3155 12 )^ 12 
+ (-8028487680 + 43726458885 2 + 9531093765 4 + 935646885 s 
+ 32762325 s + 7955 10 )^ 13 + (5278310400 - 10395422725 2 - 1832723845 4 

- 121829045 6 + 163025 8 )^ 14 - (2310758400 - 701629445 2 - 111624325 4 
+ 5909845 6 )^ 15 + (673812480 + 335319045 2 + 35189765 4 )^ 16 

- (124616704 + 90823685 2 )^ 17 + 11939840^ 18 ] 

A zV4 + 5 2 -Az + 2J(l-z)(z 2 -S 2 ) 

- 155 2 f-2 + ^ 4 ln[ V , = '} 

zy/4 + 5 2 - 4z - 2^(1 - z){z 2 - 5 2 ) 

x [230405 12 - (1920005 10 + 691205 12 )^ + (5191685 s + 5076485 10 

+ 898565 12 )^ 2 - (2027525 6 + 13378565 s + 5339525 10 + 639365 12 )^ 3 

+ (1474565 4 - 4633605 6 + 14155525 s + 2451205 10 + 268805 12 )^ 4 

- (614405 4 - 26309125 6 + 7310725 s + 131365 10 + 61445 12 )^ 5 

- (7372805 2 + 2191365 4 + 44523525 s - 1697605 s + 381445 10 - 10565 12 )z 6 



+ 
+ 



11 



+ (8355845 2 - 7249925 4 + 39598725 s - 22565 s + 103805 10 - 2525 12 )^ 7 

- (1179648 - 32604165 2 - 41876485 4 + 1 7985285 s - 460245 s - 4885 10 

- 425 12 )/ + (3244032 - 106291205 2 - 73067525 4 + 699365 s - 328045 s 

- 4655 10 )^ 9 - (3031040 - 133734405 2 - 60774405 4 - 2185925 6 - 74185 s 

- 215 10 )z 10 + (950272 - 83466245 2 - 24317445 4 - 611645 6 - 3455 8 )^ n 
+ (57344 + 24629765 2 + 4284165 4 + 35545 6 )^ 12 - (63488 + 2593285 2 

+ 222885 4 )^ 13 + (17408 - 53765 2 )^ 14 - 4096^ 15 ]}. (36) 



128tt z 2 -5 2 

455 6 s 4 (5 - z){-2 + z) w z 7 (5 + z) ^ (4 + 5 2 - 4s) (1 - z) 



+ 



x [18432005 10 (4 + 5 2 ) 2 z - 6144005 8 (4 + 5 2 ) (72 + 1015 2 + 155V 
+ 307205 6 (10368 + 428165 2 + 383685 4 + 95745 6 + 6795 8 )^ 3 

51205 6 (390528 + 8287365 2 + 5195965 4 + 1004535 6 + 54875 8 )^ 4 
51205 4 (-110592 + 10581445 2 + 15311045 4 + 7443945 s + 1151975 s + 48695 10 )^ 
12805 4 (-2955264 + 63953285 2 + 71243845 4 + 28659245 s + 3644425 s 
118895 10 )/ + 2565 2 (1 105920 - 436748805 2 + 287670085 4 + 264015445 6 
93683205 s + 10068415 10 + 249075 12 )/ - 645 2 (28139520 - 2965990405 2 
+ 616898565 4 + 448136005 6 + 158347765 s + 15165445 10 + 270595 12 )/ 
+ 
+ 
+ 

+ 



+ 
+ 



32(-2949120 + 1602764805 2 - 6171863045 4 + 567040005 6 + 101676485 s 
71124325 10 + 7181265 12 + 67535 14 )^ 9 + 16(47185920 - 5278515205 2 
7299456005 4 - 1491887365 6 + 135448325 s + 473605 10 - 1441805 12 + 24835 14 )^ 10 
8(365363200 - 10465853445 2 + 2240471045 4 - 4155874565 6 - 25762565 s 
17102005 10 + 277645 12 + 6355 14 )^ n + 4(1772748800 - 9909534725 2 

- 6652574725 4 - 6641145605 6 - 294654085 s + 2774605 10 - 211755 12 + 2105 14 )z 12 

- 8(1473757184 + 1993175045 2 - 2340190725 4 - 1502308645 6 - 75871605 s 
1243035 10 + 3005 12 )^ 13 - 4(-3491954688 - 10955100165 2 + 452657925 4 

+ 713932005 6 + 27949525 s - 185995 10 + 1055 12 )z 14 - 4(2995691520 
+ 9464064005 2 + 851501445 4 - 38959925 s + 956745 s + 4255 10 )^ 15 
+ 8(930959360 + 2437916165 2 + 243438725 4 + 7563925 s + 46175 8 )^ 16 



96(34273280 + 66155525 2 + 4429685 4 + 56375 6 ^ 17 



z 



+ 1024(956608 + 1161585 2 + 33675 4 )^ 18 - 16384(10698 + 6315 2 )^ 19 + 14581760z 20i 
+ [18432005 12 (4 + 5 2 ) - 6144005 10 (72 + 895 2 + 155 4n 



614405 8 (1296 + 42285 2 + 28195 4 + 3425 6 )z 2 
153605 8 (25440 + 428565 2 + 205935 4 + 18695 6 )^ 3 



+ 76805 6 (-25344 + 1040485 2 + 1208845 4 + 487285 s + 33915 8 )^ 4 

- 38405 6 (-283392 + 2184485 2 + 1973005 4 + 788075 6 + 42725 8 )^ 5 

+ 38405 4 (46080 - 7258885 2 + 1072645 4 + 749425 s + 444955 s + 19015 1( V 



12 



- 19205 4 (540672 - 22512645 2 + 110805 4 - 374825 6 + 353475 s + 12095 l0 )z 7 
+ 4805 2 (-98304 + 58122245 2 - 94750725 4 + 318565 6 - 3159445 8 + 404765 10 
+ 11475 12 )/ - 4805 2 (-540672 + 93859845 2 - 72807685 4 + 3622245 s 

-1793325 s + 92305 10 + 2295 12 )/ + 2405 2 (-2670592 + 201287685 2 - 88206085 4 

+ 8315685 6 - 1106965 s + 39335 10 + 855 12 )z 10 - 605 2 (-15597568 + 596500485 2 

- 184381445 4 + 17799045 6 - 848005 s + 29215 10 + 495 12 )/ 1 + 305 2 (-30146560 
+ 611061765 2 - 174074885 4 + 9806885 6 - 266365 s + 6405 10 + 75 12 )/ 2 

- 155 2 (-41222144 + 413982725 2 - 139220485 4 + 2108965 6 - 94045 s + 195 10 )^ 13 

- 155 2 (20840448 - 79544325 2 + 41658885 4 + 190405 6 + 12425 s + 75 10 )^ 14 
+ 155 2 (7962624 - 2304005 2 + 8294405 4 + 65085 6 + 695 8 )^ 15 

- 305 2 (1140736 + 1391365 2 + 484325 4 + 2195 6 )/ 6 + 2405 2 (29312 + 40645 2 
+ 3215 4 )z 17 - 19205 2 (488 + 415 2 )/ 8 + 614405V 9 ] 

zy/4 + 5 2 - Az + 2J(l-z)(-5 2 + z 2 ) 

x In — V - }. (37) 

zy/A + 5 2 -4z- 2^(1 - z)(-5 2 + z 2 ) 



1287T 



455 6 s 4 (5 - z){-2 + z) w z 7 (5 + z) \ (4 + 5 2 - 4s) (1 - z) 
x [-6144005 10 (4 + 5 2 )(4 + 35 2 )^ + 6144005 8 (96 + 1645 2 + 1035 4 + 155 6 )^ 2 



,2 _ §2 



102405 6 (10368 + 474885 2 + 394405 4 + 173025 6 + 20375 8 )/ 

+ 51205 6 (124032 + 3424005 2 + 1741485 4 + 557795 6 + 54875 s )/ 

- 51205 4 (-55296 + 2699845 2 + 7177925 4 + 2417665 6 + 571835 s + 48695 10 )z 5 

+ 12805 4 (-1480704 + 5484165 2 + 38994885 4 + 8957405 6 + 1534385 s + 118895 10 )^ 

- 2565 2 (-368640 - 213081605 2 - 103224645 4 + 181646005 6 + 28269445 s 

+ 3192835 10 + 249075 12 )/ + 645 2 (-19292160 - 1311872005 2 - 976829445 4 

+ 496544965 s + 48716245 s + 2216405 10 + 270595 12 )/ 

- 32(2949120 - 1869824005 2 - 1952512005 4 - 1891281925 6 + 583938885 s 
+ 31686245 10 - 1551985 12 + 67535 14 )/ - 16(-47185920 + 9980723205 2 

- 517463045 4 + 1328963845 6 - 737297925 s - 26174085 10 + 2697885 12 

+ 24835 14 )^ 10 + 8(-365363200 + 34124021765 2 - 8156200965 4 - 1656369925 6 

- 919752325 s - 26045525 10 + 1855165 12 + 6355 14 )^ n - 4(-1772748800 

+ 80149708805 2 - 15889187845 4 - 4626428165 6 - 812920325 s - 17141165 10 

+ 246475 12 + 2105 14 )/ 2 + 8(-1473757184 + 33237422085 2 - 3467015685 4 

- 1041226725 6 - 101993365 s - 2372655 10 + 30205 12 )^ 13 - 4(-3491954688 
+ 38653194245 2 - 632648965 4 - 367761445 6 - 25871605 s + 244255 10 

+ 1055 12 )^ 14 - 4(2995691520 - 14822348805 2 - 569283205 4 - 123445 6 

- 2443905 s + 4255 10 )^ 15 + 8(930959360 - 1503029765 2 - 76095365 4 

- 5298805 6 + 46175 8 )/ 6 - 96(34273280 + 2176005 2 + 327925 4 + 56375 6 )^ 17 



13 



+ 1024(956608 + 612945 2 + 33675 4 ),2 18 - 16384(10698 + 6315 2 )^ 19 

+ 14581760z 20 ] - [6144005 12 (4 + 35 2 ) - 6144005 10 (24 + 315 2 + 155 4 )z 

+ 614405 8 (432 + 16045 2 + 9095 4 + M25 6 )z 2 - 153605 8 (7968 + 191125 2 

+ 55355 4 + 18695 6 ),z 3 + 76805 6 (-11520 + 216485 2 + 683485 4 + 91205 6 

+ 33915 8 )z 4 - 38405 6 (-128256 - 395045 2 + 1659805 4 + 54175 6 + 42725 8 )z 5 

+ 38405 4 (39936 - 2991365 2 - 2275205 4 + 1469145 6 - 47835 s + 19015 10 )/ 

- 19205 4 (491520 - 7408645 2 - 7941205 4 + 1989705 6 - 132675 s + 12095 10 )z 7 

+ 4805 2 (-98304 + 53862405 2 - 19537925 4 - 34259045 6 + 4240085 8 - 301085 10 

+ 11475 12 )/ - 4805 2 (-442368 + 86200325 2 - 5303045 4 - 26687845 6 + 1750285 

- 93505 10 + 2295 12 )^ 9 + 2405 2 (-1753088 + 176855045 2 - 523525 4 - 32197285 6 
+ 1020325 s - 27975 10 + 855 1 2)z 10 - 605 2 (-7995392 + 480583685 2 - 15421445 4 

- 60711045 6 + 656805 s + 1355 10 + 495 12 )^ 41 + 305 2 (-11141120 + 434892805 2 

- 40924165 4 - 42556965 s + 7005 s + 6325 10 + 75 12 )^ 12 - 155 2 (-8323072 

+ 264314885 2 - 39508485 4 - 20128805 6 - 54925 s + 1735 10 )z 13 + 155 2 (196608 

+ 59842565 2 - 6735365 4 - 2728645 s - 1825 s + 75 10 )z 14 - 155 2 (2064384 

+ 14080005 2 + 931845 4 - 136205 6 + 695 8 )^ 15 + 305 2 (587776 + 1847045 2 

+ 249445 4 + 2195 6 )z 16 - 2405 2 (21632 + 41925 2 + 3215 4 )^ 17 

+ 19205 2 (440 + 415 2 )^ 18 - 614405V 9 ] 

zy/A + 5 2 -Az + 2J{l-z){-5 2 + z 2 ) 

x In — — v :}. 

zy/A + 5 2 -Az - 2yJ(l - z){-5 2 + z 2 ) 
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+ 2 flipped graphs 




Figure 1: Feynman diagrams for e + + e — > 7* — > charmonium + cc. 
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